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Abstract 

The dependence of the Casimir force on material properties is important for both futm'e appli- 
cations and to gain further insight on its fundamental aspects. Here we derive a general theory 
of the Casimir force for low-conducting compounds, or poor metals. For distances in the microm- 
eter range, a large variety of such materials is described by universal equations containing a few 
parameters: the effective plasma frequency ujp, dissipation rate 7 of the free carriers, and electric 
permittivity £00 for uj > cUp (in the infrared range). This theory can also describe inhomogeneous 
composite materials containing small regions with different conductivity. The Casimir force for me- 
chanical systems involving samples made with compounds that have a metal-insulator transition 
shows an abrupt large temperature dependence of the Casimir force within the transition region, 
where metallic and dielectric phases coexist. 

PACS numbers: ll.lO.Wx, 73.61.At 
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I. INTRODUCTION AND MOTIVATION 



The Casimir force- has demonstrated the reahty of zero-point field fluctuations, which 
played a significant role in the development of quantum field theory (see, e.g., the mono- 
graphs^i^ and review papers^i^i^i^) . The Casimir effect attracts considerable attention be- 
cause of its numerous applications in quantum field theory, atomic physics, condensed matter 
physics, gravitation and cosmology.-i^>^i^»^'^i^»^'i^ The experimental observation of the Casimir 
force is of fundamental importance. Despite the fact that the magnitude of the Casimir force 



is established by a 
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i Furthermore, this force is relevant for various 
nanomechanical devices, where the space separation of nearby plates is very smalL-i^ii^ 



A. Casimir force for good metals and dielectrics 

The Casimir force between two macroscopic samples is caused by a spatial redistribution 
of the fluctuations of the electromagnetic field compared to that of free space because of 
the presence of the samples. For the simplest case of two parallel perfectly conducting thick 
metallic plates placed in vacuum and separated by a distance /, the Casimir force per unit 
area of the sample at zero temperature can be written as 

Fc = — - (1) 
^ 240 /4 ' ^"-^ 

where c is the speed of light, and h is the Planck constant. For dielectric bodies with 
frequency-dependent dielectric permittivities, the value of this force has been found by 
Lifshitz.-i^ If the permittivity e is frequency- independent, for two equivalent dielectric bodies 
or for a dielectric sample and an ideal metal, this force can be written as 

„ tt'^ ch ( e — \ 



240 /4 



^vie). (2) 



where i/ = 2 for two equivalent dielectric bodies and v = \ for the interaction of a dielectric 
sample and a metal. The function '^yie) 1 when e ^ 1, and '^v[s) decreases when e — 1; 
in particular, ^ 1) = 0.46 and if^is ^ 1) = 0.35. Strictly speaking, equations of 

type ([I]) or ([2]) are valid when I < hc/kT, where only the zero-point fluctuations of the 



electromagnetic field are important (see Refs. 



HOP 



for details). At room temperature. 
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this inequality is valid for distances less then a few micrometers. Below we will only consider 
this range. 



B. Material aspects of the Casimir force 

To study the Casimir force, different materials can be used. Indeed, it is important to 
understand how this force is affected by the choice of different materials. For example, 
recent studies, using silicon with different degrees of doping22ii2i or materials for sensors, 
like Vanadium oxide^i, have shown numerous specific features which are absent in the good 
metals traditionally used to study the Casimir force. 

The investigation of material-dependent features of the Casimir force is important not 
only for future applications, but also for fundamental physics. To discuss the material- 
dependent aspects of the Casimir force, let us note the following. The well-known results 
present in the expressions ([1]) and ([2]) are obtained for frequency-independent values of 
the electrical permittivity e. For metals, this means e = oo for any frequency. Detailed 
investigations, taking into account the dispersion of the media, have shown^^ii^ that the 
universal formula of the type Fc oc 1//^ is valid for distances I > Aq, where Aq = c/uq, 
and ujq is the highest characteristic frequency of the media. Beyond this approximation, the 
Casimir force F can be written^'i^'iSi'i^ as 



h 



oo 
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27r2c3 



CX-$[^«)], (3) 







where e = £{iC) is the complex permittivity of the media, the summation over the Matsubara 
frequencies is replaced by integration over ( (this is adequate^^ when / < hc/kT), ^[e{i()] 
is a functional of the function e{i(), 



oo 

2, > 1 1 



$[^«)] = / ^^-^ + ^^-^ ) , (4) 
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. s + p ^ pe + s , 7. — - 

Ax = A2 = , s= ^/e + p^ -1 , 

s — p pe — s 

where x = 2p(l/c. Two terms in Eq. (jl]) describe the contributions of the modes with 
two different polarizations of the electric field, parallel to the surface and parallel to the 
incidence plane (which includes the normal to the surface and the wave vector of the photon) , 
respectively. The exponents u = 2 and v = 1 correspond to the same cases as for Eq. ([2]), 



namely, the interaction between two equivalent dispersive media (z/ = 2), and dispersive 
medium, interacting with an ideal metal {u = 1). The general properties of the function 
e{i() are the following: e{i() is a monotonic function of (, and e{() —* 1 for the values of ( 
higher than all the characteristic frequencies ( > ujq of the medium. For metals, the plasma 
frequency Up is the highest frequency Uq. Thus, the standard Casimir result ([T]) is valid for 



large distances I > \p = c/up between the plates, see Ref . [l8|. For the opposite limit case^^ 
of smaller distances, / < 



oo 



n'^0) = ^^, ^ = /(i^)".C, (5) 



where the real dispersion, e.g., the dependence of the media permittivity on the frequency, 
is used. 



C. Caviats and limitations 

It is worth noting here that, as far as we know, only one experiment-^ has been performed 
using the parallel-plate configuration originally envisioned by Casimir. Most measurements 
of the Casimir force have studied the interaction of a spherical probe with a flat substrate, 
using the so-called Proximity Force Theorem^ to relate the force for different geometries 
of the experiment to the force between two parallel plates. The experimental search for 
corrections to this approximation has been done recently.— For the original plane-parallel 
geometry, the accuracy of the measurements^ of the Casimir force, done for distances rang- 
ing from 0.5 to 6 micrometers, is not very high, within 15%. A significant difficulty has been 
the necessity to keep the samples parallel during the measurements at different distances. 
Some of these problems, in principle, could be overcome by measuring the Casimir force in 
a fixed geometry of the experiment (fixed /, for plane-parallel geometry) by varying some 
parameters of the sample. The media properties could be changed by varying the temper- 
ature of the sample. Varying the carrier density of semiconductors by laser irradiation has 
also been proposed recently.—"^ 

The Casimir force for standard metals has a weak temperature dependence. For metals, 
the Drude formula, e = 1 + ujp/({( + •y) is typically used, where Up is the metal plasma 
frequency and 7 is the relaxation rate. For typical metals like copper, aluminum or gold, 
the plasma frequency is practically temperature-independent, and the only way to modify 
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the Casimir force by changing the metal parameters is via the temperature dependence of 
7. For such metals, 7 <^ tUp, and the corresponding corrections are small. Another problem: 
for standard metals the value of \p = c/up lies in the ultraviolet region, Ap < 0.1 /im. Thus, 
to observe dispersive effects, the region I < Ap should be investigated, which is quite difficult 
experimentally.— This limitation can be overcome by using thin metallic films,— but even 
for this optimal case the temperature corrections are not higher than a few percent. 

D. Casimir force for pure metals and compounds 

Numerous compounds are known for which the carrier density and plasma frequency Up 
are abnormally small. The investigation of such conducting systems, which can be called 
"poor metals" , is of interest from the point of view of both fundamental physics and ap- 
plications. Examples include highly doped silicon^^i^, left-handed materials^^, transition 
metal oxides showing the metal-insulator transition (MIT)^, cuprate high-temperature su- 
perconductors^, and manganites where the phenomenon of colossal magnetoresistance is 
observed.— For all of these systems, both the free carrier density and the plasma frequency 
ujp are much smaller than for standard good metals. This means, that in contrast to the 
usual metals, Up is not the highest frequency of the material. The Drude behavior is ob- 
served up to infrared frequencies, but with a relatively large value of e = £00 when u 3> ujp] 
this value, £00 — 5-10, is determined by transitions of electrons in occupied bands. Thus, 
e ^ 1 within a wide frequency region, including the "metallic region" , from small u up to 
a few Up. The dissipation rate 7 for poor metals can be quite high, of the order of a few 
percent, or even a few tenths of Up. The manifestation of the dispersion for the frequencies 
corresponding to distances of the order of a few microns provides the possibility to control 
the Casimir force by varying the parameters of the metal. Recently, measurements of the 
Casimir force between a metallic sphere and a sample made with a low-conduction medium, 
like silicon with different degrees of doping and vanadium dioxide VO2, were proposed^ for 
small separations, around 200-400 nm. 

Here we develop a general theory of the Casimir force for low-conducting compounds, i.e., 
poor metals. We show that, for distances in the sub-micrometer and micrometer ranges, the 
Casimir force for a large variety of such systems can be described by formulae that depend 
on a small number of parameters, without details of the total spectral characteristics. The 
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inhomogeneous composite systems considered here, containing small regions of different 
properties, can be described within this theory. The application of these results to the 
region of the metal- insulator transition, where the metallic and dielectric phases coexist, 
produces a very pronounced temperature dependence of the Casimir force. 



II. DERIVATION OF THE CASIMIR FORCE FOR POOR METALS 



For general dispersive media, the Casimir force is determined by the integral in Eq. (jlj). 
Keeping in mind the large variety of poor-metal parameters discussed above, we now need to 
develop an analytical approach to estimate the integral (jl]) and to study the role of different 
parameters, like eoo or 7/co'p, describing the system. Let us now use a two-scale model for 
e{Q as follows: 



^(0 



1 + 



UJZ. 



(6) 



C(C + 7) 

where the function E{C,) describes the high frequency dependence of e. As we will show 
below, the detailed properties of this function are not important in the region of interest: 
/ ~ 1 yum. The function -E'(C) is almost constant, E{Q = Eoo, for all the metallic region, 
Up u Uq, and tends to one for u uq. Obviously, for such a model the standard 
Casimir behavior in Eq. ([T]) is valid at large enough distances: Z ^ Ap ~ 1 fim. 

To calculate the Casimir force for distances of the order of c/cjp ~ 1 fim we use the 
general equation ([3]) rewritten as 



h 



27r2c3 



1 + 



C(C + 7) 



00 

+ J C'dC^EiC)] 



Here the value {u) is chosen in the intermediate region: 

iUp <t: (id) ^ UJo. 

Therefore, we replaced E{() by Eoo in Eq- © for the first integral and omitted the Drude 
multiplier for the second integral. 
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Expanding the integration region over ( in both integrals up to < C < oo, and sub- 
tracting the extra terms, we present the Casimir force in the form. 



F = + AF, 



with 



pirn) 



h 



27r2c3 



AF 



h 



27r2c3 



(7) 



(8) 



(9) 



In the frequency region, u c/l, which is an important regime for the functions 

$[£'((■)] and $[£oo] in Eq- © almost cancel each other. Therefore, the term AF is relatively 
small. A more detailed analysis gives 

AF = — ^ < F. 

Thus, in the region of interest, Aq ^ / ~ Xp, the Casimir force is described by the first term 
in Eq. (^i), F = F^"^). 

Now we introduce the variable z = (l/c and write the main contribution to the Casimir 
force as 

(10) 



F = -i?.n^Fc-n 



240 

where Fq is the Casimir force ([T]) for ideal metals, the prefactor 11 depends only on the 
dimensionless parameters / = l/Xp, Eoo, and a = ■y/ujp, 

oo oo 

^ 120 /■ o , /■ 2 , 
n = — - I z^dz / p dp 
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A^exp(x) — 1 A2exp(a;) — 1 
where Ai and A2 are given by Eq. (j4]) with 



(11) 



e = Er 



1 + 



z{z + al) 
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FIG. 1: The normalized Casimir force 11 = F/ Fq versus l/\p for some values of parameters (shown 
near curves). This plot compares the results of numerical calculations of the two-dimensional 
integral Eq. (jlip for the prefactor IT (symbols) with the results within the approximate approach 
based on Eq. (I12p (solid curves). This plot shows a very good agreement between both. 

A. Computing the Casimir force integrals 

To proceed further, let us change the variable z hy x = z/2p m Eq. ffTTl) . and note, that 

oo 

the integral / [Ae^ — l)^^x'^dx, with A = A{x,p) being a smooth function of x, can be 



approximated by l/A{xo,p) with xq ~ 4. The problem is then reduced to calculating two 

one-dimensional integrals, Ji and J2: 

00 

1 

where Ai and A2 are given by Eq. (4) with the substitution 

(13) 

The validity of this approximation is confirmed by the numerical calculation of integral (|TT1) . 
as shown in Fig. 1. The function n(/) found numerically is shown in Fig. [21 A simple analysis 
of Eq. (fT2l) gives us two limit cases. 

For small / < 0.3 the value of a plays a minor role. In this region, the value of 11 
practically does not depend on I and reproduces well the Lifshitz's result ([2]) for dielectric 
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FIG. 2: (Color online) The normalized Casimir force 11 = F/Fq versus the parameter I = l/Xp, 
for v = 1,2, using the typical value £oo = 8, as well as the smaller Eoo = 1-5, and different values 
of the dissipation parameter a = 'j/iOp. The horizontal line on top gives the asymptotic value 
F/Fq = n = 1 for an ideal metal. 

media with a C- independent e = Soo and 7 = 0, 



We now emphasize that the dependence of the Casimir force, proportional to Co/l^, see 
Eq. ([5]), is not realized for any Eoo ^ ^■ 

Otherwise, in the hmit case / — oo, the integrals Ji = J2 = 1, and the ideal Casimir 
limit ([T]) is recovered. In contrast to the case of small values of /, the dependence of 11 on / 
for large, but finite values of / shows an interesting and unexpected behavior: the approach 
to saturation is quite slow, especially for large values of 



In other words, it is hard to reach the metallic limit value of 11 = 1 when a > 0.1, for the 
most interesting region I < 10. 

To understand this behavior, let us now investigate in more details the factor 11 for not 
so small values of /. As has been mentioned above, it is a sum of two contributions from 
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FIG. 3: Integrals Ji and J2, defined by Eq. (jl2p . describing the Casimir force versus l/\p for 
£00 = 8 and different values of tlie dissipation parameter a = "y/ujp (shown in the figure). Symbols 
depict the results of numerical calculations. 

the electromagnetic fields of different polarizations. It is convenient to examine the first 
and second integrals separately. Numerical calculations show that the behavior of these two 
integrals, Ji and J2, is essentially different for the same values of parameters, as shown in 
Fig. El 

The two interesting features (i.e., the slow approach to saturation and the essential de- 
pendence of n on a = '~i /ojp) are mostly associated with the first integral, Ji, which describes 
the contribution of the fiuctuations with the electric field parallel to the surfaces of plane- 
parallel samples. This integral Ji can be calculated analytically. For a = 0, it can be written 

as 




H , ■ arctan \ , 14 

6Va^^ y a + 1 v/^2:^ + 5+ I ' V y 

where we introduce the notation 

P 

a^ = l + £oo^, 6^ = £00 - 1- 
For non-zero dissipation rate, a 7^ 0, the analytical formula for Ji is very long and incon- 
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FIG. 4: The integral Ji versus l/\p for = 8 and different values of the dissipation parameter 
a = 7/wp. Symbols depict the results of numerical calculations. The detailed behavior of Ji for 
small / is present in the inset. The dotted line describes the analytical result Eq. (|14p for a = 0, 
solid lines are drawn in accordance with the approximate formula Eq. (jlSp . 



venient for real estimates. But for the case of interest, / > 3, the integral Ji can be well 
approximated with the simpler expression 



appr 



(3 + £oo)(3 + SaO + ^ej"^ - 2\/3 + 5a[ 
(3 + £oo)(3 + SaO + 3^00^^^ + 2\/3 + bal 



(15) 



as shown in Fig. H] 

This equation explains the complicated behavior of the factor FI = F/Fq for large /, and 
especially, the role of the dissipation constant a. For any finite value of a and extremely 
large /, / 1 and / 1/a, Ji versus / has a very slow inverse-square-root dependence, 

AV5a 1 



Jl - 1 



/ > 1, 1/a. 



(16) 



For very small a -C 1, the intermediate region 1/a ^ / ^ 1 appears. For this region, the 
behavior of Ji is sharper. 



2 1 , r 

Jl ~ 1 — ■ - , 1/a > / > 1 . 



(17) 
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FIG. 5: The integral J2 calculated numerically for Soo = 8 and £00 = 1-5 and different values of a 
(symbols). The shapes of the curves depend very weakly on a. 

An analytical expression for J2 in terms of elementary functions cannot be written. For- 
tunately, the integral J2 exhibits a simpler behavior than Ji, and for its description we can 
use a simple approximation. First note the quite weak dependence of the shape of the func- 
tion J2(0 ^he value of a, as shown in Fig. [5l For the regime of interest here, Eoo ^ 1, the 
difference between the values of J2 for a = 0.3 and a = is maximum near the range I ~ 
(10-15), and does not exceed 3%. Indeed, all the curves with < « < 0.3 merge together, 
and for describing of J2 within an accuracy of 1.5%, the function ^2(0 found for a = 0.1 
can be used. Even for small £00 = 1-5, the inaccuracy of this approximation is less than 5%. 

Numerical data are well fitted by the very simple formula 

J„, = . (18) 

where t] ~ (0.5-0.6), J2,l determines the value of J2 for small / <^ 1, as shown in Fig. [61 
The quantity J2,l describes the contribution of J2 to the Lifshitz's result ([2]) for dielectric 
media with e = Eoo and a = 0. 

Thus, we can present a simple description of the second integral J2: it is practically 
independent on the dissipation parameter a, and the dependence on ^oo is governed only by 
the Lifshitz contribution J2,l- The asymptotic behavior of J2 at large distances / ^ 1 is of 
the form 1 — (1 — J2,l)/'7^5 which is much weaker than the inverse-square-root dependence ( 1T6|) 
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FIG. 6: The integral J2 calculated numerically for a = 0.1 and very different values of the dielectric 
permittivity, £00 = 8 and £00 = 1-5 (symbols). The curves describe the approximating function 
Eq. (jlSp . Note the very good agreement between numerical data and approximating functions. 

for Ji. For large Eoo ^ 1, when (1 — J2,l) oc l/^oo ^ 1, the dependence ^2(0 is especially 
weak, even compared with that for Ji in the intermediate region (JTTI) . 

B. Change in the Casimir force near the metal-insulator transition 

The analytical formulae derived above give a good description of the behavior of the 
Casimir force when the metal-insulator transition occurs. Usually, the metal-insulator tran- 
sition is associated with an abrupt change, by a few orders of magnitude, of the conductivity 
at a transition temperature T = Tc. Let us start with a rough picture, assuming that a 
metallic phase has a finite value of the plasma frequency whereas for the dielectric phase 
the plasma frequency is zero. Using the results obtained above, one can expect a drastic 
change of the Casimir force between two plane-parallel samples caused by the change of the 
parameter Ap, very near the metal-insulator transition. 

We stress that the change of the force is not connected with changing the physical distance 
I, but with changing the dimensionless quantity I = l/Xp, caused by the change of the plasma 
wavelength Ap = c/up. Thus, one can expect a jump-like behavior of the Casimir force when 
changing the temperature across Tc. Within the transition region, the force changes from 
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the "metallic" value , typical for finite values of /, to the very different value F^, for an 
insulator when / <^ 1. 

The important quantity here is the change of the Casimir force, AF = {F^ — F^). To 
estimate F^, we can use the Lifshitz formula ([2]) valid in the limit / <^ Xp, which corresponds 
to the dielectric phase. The value of F-^ in the metallic phase corresponds to large, but finite 
values of l/Xp. To estimate note that the dependence of 11 on / at 10 > //Ap > 2.5 is 
mainly provided by the integral Ji, whereas J2 can be replaced by one. Thus, the concrete 
value of the coefficient t] in Eq. f|T8l) is not important. Combining all these data together, 
and restoring the initial parameters of the media, Up and 7, we arrive at the simple estimate, 

vr^ ch 



AF 



240 



1 - 



1 



V^i(£^oo) 



2 



5c7 



(19) 



- 00 I -^7 

where the function ipi{e) describes the Lifshitz's result for the interaction of a dielectric 
sample and an ideal metal. When the value of 7 is small enough, as for manganites, for the 
distance / such that / <^ c/7, Eq. ( |T71) is valid, and the formula for AF reads 

TT^ ch 



AF 



240 



1 



1 



+ 1 



(20) 



21 



Note that our results differ significantly from the theoretical estimates given in Ref. 
In particular, the value of AF in Ref. 2ll is proportional to the temperature T. The linear 
dependence AF on T can be expected for large enough separations / ^ ch/kT, that is, 
larger than a few microns, and cannot appear for small separations. 

It is worth noting that the relative change of the force 

F> - F< 
AF = 

F> 

is larger for long distances /, when the value F^ of the force for media in the conducting 
state is larger than the limit value F^ describing the case of small ujp and small l/Xp. This 
feature is determined by the quite slow change of the function n(/) at not so small values 
of £00, as shown in Fig. [2l 



III. COMPOSITE MEDIA AND THE INTERMEDIATE REGION FOR METAL- 
INSULATOR TRANSITION 

The very abrupt (by a few orders of magnitude) change of the conductivity at the metal- 
insulator transition occurs for the dc case only. At finite frequencies, the behavior of the 
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FIG. 7: (Color online) The effective permittivity ecfT(^C) (iii units of £00)5 for d = 3 and different 
concentrations / of the metal phase, as a function of the dimensionless variable X"^ = (^(C + 7) /ujp. 
Inset: effective plasma frequency cOp^cS (in units of ujp) versus /, in the coexistence region. 

complex permittivity of compounds near metal-insulator transition is more complicated and 
the jump-like behavior, typical for the static conductivity, does not arise for e = e{i(). 
Within the finite transition region, the presence of a non-uniform state with coexisting 
metallic and insulator phases is well established for all systems showing a metal-insulator 
transition. Obviously, this effect is of great interest for studying the Casimir force. The 
effective-medium approach suggests that the metallic and insulating regions coexist as in- 
terpenetrating clusters, providing a percolation picture^^ of the metal-insulator transition at 
T = Tf.. When the transition is of first order, the phase-separated regions are mesoscopic, in 
the 100 nanometer range, and quasistatic objects (giant clusters) have approximately equal 
electron densities. 

To describe the Casimir force for such a nonuniform state, we have used the effective- 
medium approximation^™ developed for composite metal-insulator media. This approxima- 
tion has been used for explaining the optical properties of VO2 near the metal-insulator 
transition.— In this model, the effective value of e = EesiC) is determined by the concentra- 
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FIG. 8: Casimir force for the interaction between two equivalent poor metals {u = 2) (with s^o = 8 
and a = 0.3) versus the concentration / of the metallic phase in the coexistence region. The force 
is normalized by values of the force Fq for ideal metals in the same geometry, and for different 
values of l/Xp, where Xp is determined by the plasma frequency in the pure metallic phase. 



tion / (0 < / < 1) of the metal phase following the equation, 

^eff , /-I r\ ^eff 



/ 



+ (1-/) 



0, 



(21) 



em + {d- l)ecs £i + {d- l)ees 

where Em and Ei are the frequency-dependent permittivities for the metallic and insulating 
phases, respectively. Also, d — 2 and d — Sior the thin film (thickness smaller then the grain 
size) and bulk sample, respectively. In the intermediate region, the effective permittivity 
EesiK) as a function of the phase concentration / can be written as follows: 



2^^ = (2/-l) 



and 




(22) 



C(C + 7) 



(23) 



6(l + /)a;2 



(3/ - iH 

C(C + 7) 



for 0? = 3 . 
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These equations predict an infinite value of £eff(^C) wlien C (tliat corresponds to a 
metallic conductivity) when /c < / < 1 only, where /c = l/rf is a percolation threshold, see 
Fig. El Otherwise, a dielectric behavior is present, with a finite value of £efr(^C) when C 0, 

SesiC = 0) = Y^fd ^ ^ d ' 

as shown in Fig. [71 In the metallic region (above the percolation threshold, for / > /c), the 
behavior of ecs{iC) at small ( is determined by the effective plasma frequency ujp^es, 

Ees — > goo . , N when C — ^ 0. 
C(C + 7) 

The value of cu^ increases linearly with / from zero, at / = /c, until cj^, at / = 1. Thus, 
a square root behavior of the effective plasma frequency ci^p efr over (/ — /c) is present in the 
metallic region, see inset in Fig. [71 

It is useful to note here that a linear temperature dependence of uJ^ was observed^^ in 
Lao.ySro.aMnOs for T < T^. Thus, we can describe the Casimir force in the intermediate 
region as a series of curves, with their shape only depending on //Ap, as shown in Fig. [HI 



IV. PREDICTIONS FOR SPECIFIC MATERIALS 



Using the results obtained in the previous sections, here we estimate numbers for different 
materials showing the metal-insulator transition. To study the Casimir force in the vicinity 
of the metal-insulator transition, we choose two typical compounds: vanadium dioxide VO2 
and the manganites exhibiting colossal magnetoresistance. For these two materials, the 
general Drude behavior of permittivity, with typical values of Ap of the order of 1 /im and 
with relatively large values of e, £00 ~ 5-10, is observed in the infrared region of interest. 



A. Vanadium dioxide VO2 

Vanadium dioxide, VO2, shows a jump in the static conductivity (a metal-insulator tran- 
sition) a little bit above room temperature, at T = Tc ~ 68 °C. The pure metallic phase 
of VO2 is realized at T > 88 °C, and pure insulator phase^^ (more exactly, semiconducting 
phase with a gap of the order of 1 eV) at T < 60 °C. For vanadium dioxide, the phase sepa- 
rated state has been observed^^ within a finite temperature range, between 60 °C and 88 °C, 
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by measuring the optical properties of VO2. Recently, such state was directly observed^ 
via scanning tunneling spectroscopy. For all temperatures where the metallic conductivity 
is present, the generalized Drude behavior is observed up to the infrared frequency, with a 
relatively large value of ^oo — 9 and Xp = c/ ujp = 1 fim. The phonon contribution to the 
value of e, typical for the infrared region, is screened by free carriers, and the value of £00 — 9 
is kept until the high-frequency region, with A > Aq ~ 0.1 fim, where the value of {e — 1) 
vanishes. The value of the dissipation rate 7 for this compound is large enough, 'y/up 0.3 
for VO2, and the data for large a should be considered. 

For VO2, the Casimir force increases when increasing the temperature through the tran- 
sition region, from 60 °C until 88 °C, see Fig. IHl^a). The value of a is quite high, and the 
calculated change of the Casimir force is essentially smaller than for the naive estimate as 
the difference between the values Fq for an ideal metal and Fl for a dielectric, see Eqs. ([1]) 
and ([2]). The relative change of the Casimir force is maximal for large enough distances, 
e.g., / ~ 4Ap ~ /xm. 

B. Manganites. 

Manganites (with antiferromagnetic insulators LaMnOs or NdMnOa as parent com- 
pounds, after substitution of La by divalent ion) show a metal-insulator transition at the 
dopant concentration x ~ 0.3, with a ferromagnetic metallic phase in the low temperature 
range.— These systems are very popular now in the context of colossal magnetoresistance, 
based on the possibility of the metal-insulator transition induced by an external magnetic 
field, that is caused by the ferromagnetic ordering of the metallic phase. On the other hand, 
the standard temperature-induced metal-insulator transition is possible for such materials 
as well. For example, the typical compound Lao.7Cao.3Mn03 demonstrates a metal-insulator 
transition at T = Tc = 250 K. The phase separation state is present for all temperatures 
below the transition point, and the typical linear dependence of i^p^es been observed^ 
in this region. Note that this metal-insulator transition is accompanied by ferromagnetic 
ordering. In principle, it could produce an extra-force of magnetic origin near the transition 
(antiferromagnetic ordering, present for some metal-insulator transition, does not produce 
any source of long-ranged interactions). However, for large enough plane-parallel samples, 
the magnetic flux lines are closed inside the magnetic sample, and should not produce any 
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FIG. 9: Predicted temperature dependence of the Casimir force for the interaction between a 
poor metal and an ideal metal, calculated for VO2 (a) and for Lao.rCao.sMnOs (b). The force 
is normalized by its value Fc for an ideal metal, for different values of l/Xp=l, 3, 5, where Xp 
is determined by the plasma frequency in the pure metallic phase. The horizontal dashed lines 
indicate the limit values of the force for the pure insulating phase and (for VO2 only) its metallic 
phase. The corresponding values of Fc are described by Eq. ([1]); for Xp ~ 1.2 ^m these are 
Fc = 0.6 • 10-3, 0.74 • 10-^ and 0.965 • 10"^ Dyn/cm^ for l/Xp=l, 3, 5, respectively. 

serious parasitic effects. For these compounds, Up is small and the corresponding Ap ~ 1 /im. 
The main specific feature important for us here is the low value of the dissipation rate: typ- 
ical values of 'j/ujp are ~ 0.02-0.05, and the low-7 behavior of the curves shown in Figs. [21 3 
are adequate. 

For Lao.rCao.sMnOs, the metallic phase corresponds to the low-temperature range, and 
the value of the force increases when decreasing the temperature, which leads to an opposite 
temperature behavior of the Casimir force, compared to VO2. The value of a for this 
compound is relatively low, and the temperature dependence of the Casimir force is sharper 
than for the previous example. One more specific feature is the presence of phase separation 
in the whole region of the metallic phase existence. Thus, one can expect an essential 
dependence of the Casimir force for all temperatures below the transition temperature, see 
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Fig.EKb). 



V. CONCLUSIONS 

The Casimir force depends on the materials used, and we have studied some of these 
material-dependent aspects. The Casimir force Fq for a mechanical system containing com- 
pounds with a metal-insulator transition shows an abrupt temperature dependence in the 
transition region. The relative change of the force AFc when crossing the transition region 
can be quite large, of the order of the force itself for a distance ~ 5-6 microns. The relative 
change AFc of the Casimir force is larger for large distances, where the absolute value of the 
force is small, but even for a distance / = Ap = 1.2 yum it reaches 30%. The dependence of 
the force on temperature is sharp near the percolation threshold, where the static metallic 
conductivity appears. 

When measuring such tiny forces, the exclusion of any parasitic effects, like electrostatic 
forces, is essential. To avoid electrostatic forces, the usual highly-conducting samples are 
short-circuited.— This method might appear to be ineffective for the metal-insulator tran- 
sition compounds near the insulating region. However, such compounds are more semicon- 
ducting than insulating in this region and the conductivity is non-zero at room temperatures. 
Thus, we believe that the same technique could be used. To increase the conductivity in 
the semiconducting region, the usual doping by donor or acceptor impurities could be used. 
Finally, we note that the metal-insulator transition is sometimes accompanied by structural 
phase transitions, which could lead to some lattice distortions. Thus, care should be taken 
to choose materials and operating conditions that avoid these additional difficulties. 

For measurement of the Casimir force for samples made with usual metals, small separa- 
tions are preferable. The creation of experimental set-ups with very small (sub-micrometer) 
distances between samples is a serious challenge for experimentalists. As follows from our 
analysis, distances / larger than the plasma wavelength Ap are preferable for the experimen- 
tal observation of the effects, we predict around the region of the metal-insulator transition. 
For the compounds discussed above, this means distances of the order of (2-4) /im. In the 
planned experiments^^ for measuring the Casimir force using Vanadium oxide samples, the 
separations are (0.2-0.4) /im, which equals (0.1-0.3) Ap. These values are much smaller than 
the optimal values noted above. For separations / of the order of (0.1-0.3) Ap, the Casimir 
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force should follow low-/ asymptotics for any temperature (in both the metallic and the 
insulating phases). The temperature dependence of the Casimir force should be weak and 
the manifestation of the metal-insulator transition should be minor for such an experimental 
set-up. 
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